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Abstract 

We present a systematic study of the atomic and electronic structure of the Si(lll)-(5x2)-Au 
reconstruction using first-principles electronic structure calculations based on the density functional 
theory. We analyze the structural models proposed by Marks and Plass [Phys. Rev. Lett. 75, 
2172 (1995)], those proposed recently by Erwin [Phys. Rev. Lett. 91, 206101 (2003)], and a 
completely new structure that was found during our structural optimizations. We study in detail 
the energetics and the structural and electronic properties of the different models. For the two 
most stable models, we also calculate the change in the surface energy as a function of the content 
of silicon adatoms for a realistic range of concentrations. Our new model is the energetically most 
favorable in the range of low adatom concentrations, while Erwin's "5x2" model becomes favorable 
for larger adatom concentrations. The crossing between the surface energies of both structures is 
found close to 1/2 adatoms per 5x2 unit cell, i.e. near the maximum adatom coverage observed 
in the experiments. Both models, the new structure and Erwin's "5x2" model, seem to provide 
a good description of many of the available experimental data, particularly of the angle-resolved 
photoemission measurements. 

PACS numbers: 73.20.At, 71.10.Pm, 79.60.Jv, 81.07.Vb 
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I. INTRODUCTION 



The low-enerffv electronic spectrum of a one-dimensional metal is dominated by collective 
a„d c,...e excitation. Th. i. . co„.a. wit. t.e be,.avio. ot t.p.a, .et- 

als, that can be understood in terms of independent particle-like excitations usually called 
quasiparticles. These predictions are clear and well established. However, the observation 
of the Luttinger liquid behavior in real systems has proven quite elusive. One of the reasons 
for this might be that one- dimensional metals are in principle unstable with respect to the 
Peierls distortion that drives them into an insulating ground state ^. A possible route to 
avoid this limitation is the fabrication of metallic chains absorbed on surfaces. The hope 
is that the rigidity of the substrate will make the energy cost for the structural distortions 
too large and, therefore, the one-dimensional chains would remain metallic. Semiconductor 
surfaces are specially attractive for this purpose. The existence of an energy gap prevents 
the coupling of the electronic states of the chain in the vicinity of the Fermi level with the 
substrate. The one-dimensional character of these states is thus preserved. It is in this 
context, that the fabrication of monatomic wires of metal atoms on silicon substrates has 
attracted much attention in recent years. 

Monatomic wires of gold atoms are spontaneously formed on fiat and vicinal Si(lll) 
surfaces after the deposition of gold in the sub-monolayer regime (see Ref. s] and references 
therein). Of particular interest are the vicinal substrates, where gold wires run parallel to 
each step-edge and the coupling between the chains, and thus the one-dimensional character 
of the electronic states, can in principle be controlled by changing the miscut angle. Ex- 
amples of these systems that have attracted much attention in the last few years are the 

Si(557)-Au mm 
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and the Si(553)-Au [s, 3,[lJ,[l^ reconstructions. The pho- 
toemission spectra of these surfaces close to the Fermi energy are dominated by bands with 
a strong one-dimensional character that exhibit several interesting phenomena including pe- 
culiar splittings associated with the spin-orbit interaction j9| 
signature of spin-charge separation j^), fractional fillings [1 



although first interpreted as 
and metal-insulator transi- 



tions 



IjJ. The analogous to these one-dimensional structures in the case of the fiat Si(lll) 
surface is the so-called 5x2 reconstruction. 

The deposition of gold in the monolayer (ML) range on the fiat Si(lll) surface results 
in a variety of phases [16], such as x -\/3R30°, 1x1 and 5x2. The 5x2 phase occurs 
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at ~0.4 ML gold coverage 



0. 



It was first discovered about thirty years ago 



00,0 



and has been investigated using many experimental techniques since then. This includes 



low energy electron diffraction (LEED) studies 
standing wave analysis 



3,0 



20|, x-ray diffraction 2l| and x-ray 



photoemission spectroscopy (ARPES) 



22| . scanning tunnelin g ra icroscopy (STM) j2c 



26, 



27, 



25| , angle resolved 



29|,|30| and inverse photoemission 



31 




and 



high resolution electron microscopy (HREM) combined with heavy-atom holography 

Already the first structural models, based on LEED measurements, considered two atomic 
gold chains per 5x2 unit cell running in parallel 3 0|- This was later confirmed by 
HREM j^^l and seems to be firmly established. The gold chains run along the [110] and 
equivalent directions (parallel to the x2 periodicity of the unit cell). Therefore, three dif- 
ferent domains are possible for the 5x2 reconstruction on the flat Si(lll) surface. Single- 
domain surfaces, necessary for ARPES, can be fabricated using vicinal surfaces with a slight 
cut-off angle P. I23I. The presence of one-dimensional structures in this reconstruction has 
also been confirmed by the ARPES studies. Early studies found a strong anisotropic sig- 
nal near the Fermi level 3, 0] , but no evidence of Fermi- level crossing for this band was 
found 0|. Later studies at low temperature found a one-dimensional band with a strong 
dispersion along the direction of the gold chains 0, 28 1. The top of this band appears 
near the 5x2 zone boundary and disperses downward, reaching its minimum close the 5x1 
zone boundary. This band has been reported to change its dimensionality from strongly 
one-dimensional near the Fermi energy to two-dimensional at lower energies [0. In these 
studies a gap of ~0.3 eV was also identified for this band. The presence of this gap and 



its apparent closing with increasing temperature was related to a Peierls instability 



More recent ARPES results 



29 



I, both at low and room temperatures, have been able to 
identify some additional surface bands. However, the metallic or semiconducting character 
of the surface is still a matter of debate. In fact, it has been proposed that the metallic or 
semiconducting character can depend on the concentration of silicon adatoms 0,0, and 
even that semiconducting and metallic segments can alternate along the gold chains in the 
surface 0]. 



The STM images are characterized by the presence of bright, irregular protrusions |24l. |25| . 



and "Y" -shaped features |23|, |35[ with a well defined orientation respect to the underlying 
substrate. The protrusions have been established to be silicon adatoms 36], which are 
present on the surface with an optimum coverage close to 1/4 adatoms per 5x2 unit cell. 
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In spite of all these experimental studies, the structure of the Si(lll)-(5x2)-Au recon- 
struction has not been completely established yet. Earlier structural models only considered 
the adsorption sites of the gold atoms. Many of them could be ruled out on the bases of 
more detailed STM studie s [24 1 a nd the knowledge of the exact gold coverage fl^]. A few 
more refined models exist [25|, |32| . They consider both the position of the gold atoms on 
the substrate and the rebonding of the silicon atoms in the surface layer. Probably the most 
detailed structural model proposed to date is the one by Marks and Plass (MP) The 
MP model is based on a combination off-zone HREM, transmission electron diffraction and 
heavy-atom holography data. 

The first theoretical studies using first-principles electronic structure calculations ap- 
peared only quite recently. This is due to the complicated structure and the large unit 
cell of the Si(lll)-(5x2)-Au reconstruction. Kang and Lee 37| studied the MP and the 
Hasegawa- Hosaka-Hosoki (HHH) Q models using density functional theory. Their main 
conclusion is that both models fail to reproduce some of the key features of the STM images 

n 

and the experimental band structures. Using a similar methodology, Erwin proposed and 
studied new structures which are characterized by the presence of the so-called honeycomb- 
chain silicon structure jssl. One of these models (the so-called "5x2" model) seems to fulfill 
many of the constraints imposed by the empirical evidence. An interesting point raised by 
Erwin is that of the crucial role played by silicon adatoms in the stabilization of the different 
structures. According to Ref. 123, the surface energy of Erwin's "5x2" model is minimized 
for an optimum adatom coverage in agreement with recent experimental reports For 
lower adatom coverages other structures compete in stability. This is a very interesting 
result which, however, is based on approximate calculations. Due to the large size of the 
supercells necessary to simulate explicitly the effect of the different adatom concentrations, 
Erwin assumed that the main role played by the adatom is to dope the surface with elec- 
trons. He then analyzed the behavior of the total energy as a function of the number of 
extra-electrons in the substrate, obtaining a minimum for ~0.25 electrons per 5x2 unit cell. 

In this work, we present a comprehensive study of the atomic and electronic structure 
of different models of the Si(lll)-(5x2)-Au reconstruction using electronic structure calcu- 
lations based on the density functional theory. We have used two different methodologies. 



the S 
code 



ESTA code 



4G 



41[ using a basis set of localized atomic orbitals and the VASP 



42L I43| using a basis set of plane-waves. We analyze the MP model j^, the models 
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proposed by Erwin |33|, and a new model that we found during our structural optimiza- 
tions. We study in detail the energetics and the structural and electronic properties of the 
different models. We also calculate the change in the surface energy as a function of the 
content of silicon adatoms for the two most stable models. In order to do so, we perform cal- 
culations for large supercells containing realistic concentrations of adatoms: 5x4, 5x6, and 
5x8 supercells. Our new model is the most favorable in the range of low adatom concentra- 
tions, while Erwin's "5x2" model becomes favorable for larger adatom concentrations. The 
crossing between the surface energy of both structures occurs close to 1/2 adatoms per 5x2 
unit cell, i.e. near the maximum adatom concentration observed in the experiments. Both 
models, our new structure and Erwin's "5x2" model, seem to provide a good description 
of most of the experimental data. Particularly, we find a general agreement between the 
calculated and measured band structures along the direction parallel to the gold chains. 



II. CALCULATION METHOD 



39 



40 



code. We have 



45 



46| (the 



Most of our calculations have been performed using the SIESTA 
used here the local approximation (LDA) to the density-functional theory 
generalized gradient approximation (GGA) has been also used for a few test calculations), 
Troullier-Martins pseudopotentials|42(], and a basis set of numerical atomic orbitals obtained 
from the solution of the atomic pseudopotential at slightly excited energies 0, 13, Q ■ 
We have used an energy shift j40|,|49| of 200 meV. The corresponding radii are 5.3, 6.4 and 
6.4 a.u. for the s, p and d orbitals of Si, 6.2, 6.2 and 4.5 in the case Au, and 5.1 for the 
s and p states of H. The bidimensional Brillouin-zone (BZ) samplin g |50l| contained 4x4 
points for the 5x2 unit cell (and a consistent sampling for other cells j6l|) and the fineness 
of the real-space grid used to compute the Hartree and exchange-correlation contributions 
to the total energy and Hamiltonian matrix elements was equivalent to a 100 Ry plane-wave 
cutoff. This guarantees the convergence of the total energy, for a given basis set, within 
~20 meV/Au (~0.5 meV/A^). 

We modeled the surface using a finite slab, similar to that depicted in Fig. ^ For 
most calculations the slabs contained three silicon bilayers (the one at the surface and two 
underlying silicon bilayers) plus an additional layer of hydrogen atoms to saturate the silicon 
atoms in the bottom of the slab. This removes the surface bands associated to the bottom 
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surface from the energy-range of interest, i.e. from the band-gap region. We have checked 
the convergence of the results using thicker slabs for the most stable structural models of 
the surface. We use periodic boundary conditions in all three directions. A vacuum region 
of 15 A ensures negligible interactions between neighboring slabs. During the structural 
relaxations the positions of the silicon atoms in the bottom layer were kept at the bulk ideal 
positions. Unless otherwise stated all other degrees of freedom were optimized until all the 
components of the residual forces were smaller than 0.04 eV/A. To avoid artificial stresses 
the lateral lattice parameter was adjusted to the theoretical bulk value calculated using 
similar approximations to those utilized in the slab calculations, i.e. the same basis set and 
grid cutoff, and a consistent k-sampling. The values are, respectively, 5.48 and 5.42 A with 
double-C (DZ) and double-^ polarized (DZP) basis, to be compared with the experimental 
value of 5.43 A. 

Due to the large number of atoms (~70 atoms for typical slabs and up to 273 for the 
largest ones) and to the need to perform geometrical optimizations for many different struc- 
tural models, we have decided to use a DZ basis set for silicon in most of our calculations. 
This basis set includes two different functions (i.e. two different radial shapes) to represent 
the 3s orbitals of Si, and another two for the 3p shell. We have tested the performance of 
this basis set and the other parameters of our calculations using the well known Si( 11 1)2x2 
adatom reconstruction as a test. Our results using a DZ basis and a slab containing three 
bilayers reproduced very well the results for the geometry and energetics given by previous 



calculations 



5l|. In fact, using a more complete DZP basis set, which includes a shell of d or- 



bitals, or increasing the thickness of the slab (up to five bilayers) did not change appreciably 
the results. 

A scalar-relativistic pseudopotential similar to that utilized in Ref. |52 have been used for 



gold. The gold basis set included double and polarized 6s orbitals (i.e. two different radial 
shapes to describe 6s orbital plus a 6p shell) and a single 5d shell. We refer to this basis set 
as DZPs-SZd. This basis set was already used in calculations for several gold clusters, where 
it was shown to lead to results in very good agreement with more complete basis sets 5^. 
We have also tested that these basis set and pseudopotential yield to the correct bulk lattice 
parameter and band structure. 

In this work we study the relaxed structures and the energetics of several models of the 
Si(lll)-(5x2)-Au surface reconstruction. The energy differences between different mod- 
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els are of key importance since we would like to determine the most plausible structures. 
Whenever it is necessary to compare the energies of structures containing different numbers 
of silicon atoms, the silicon chemical potential is set to the total energy of bulk silicon at the 
equilibrium lattice parameter. This choice is justified by the fact that the surface should be 
in equilibrium with the bulk. A summary of our results can be found in Table |1J One can 
see that the relative energies are quite small in some cases. However, they are larger than 
the estimated error bar for the total energy (see above). Furthermore, the relative energies 
usually exhibit a faster convergence than the total energy of a single structure. It is also 
necessary to check the convergence of the results as a function of the slab thickness and the 
completeness of the basis set. Table ITTl shows the results of these tests for the most stable 
structures. In one case, the slab thickness was increased by one silicon bilayer while, in 
the other, a DZP basis set was used for the silicon atoms. In both cases the systems were 
relaxed. The results are quite stable against the change of the slab thickness. In particular, 
the energy order of the structures is not changed and the variation of the relative surface 
energies is smaller than ~0.5 meV/A^ in all the cases. The variations with the size of the 
basis set are somewhat larger. From the results in Table |n] we can estimate an error bar 
smaller than 2 meV/A^ for the relative surface energies of the different structures calculated 
using SIESTA. 

In order to check the accuracy of our predictions we decided to perform calculations for 
some of the systems with another electronic structure code that utilizes a different methodol- 
ogy. We used the VASP code j3, Q] for this purpose. We used projected- augmented- wave 
potentials and a well converged plane- wave basis set with a cutoff of 312 eV. All struc- 
tures were relaxed (the equilibrium lattice parameter of bulk silicon obtained with VASP is 
5.41 A). In Table HH we can see some of the results obtained with VASP. They are in good 
agreement with the SIESTA results, especially with those obtained with the more complete 
DZP basis set. The order between our more stable models is preserved, although the en- 
ergy difference is somewhat decreased. In particular, the new structural model found in 
the present work (model N in Tables H] and |n} is confirmed to be the most stable surface 
reconstruction between those studied here. It is also interesting to note that the energy as- 
sociated with the addition and removal of adatoms for a particular structural model seems 
to be quite independent of the details of the calculation. 

The surface BZs of the studied systems are shown in Fig. |2](a). For the 5x1 system the 
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BZ is a stretched hexagon while, for the remaining periodicities, the hexagons are distorted. 
We plot the electronic band structures of the different models along the F- ZBx2-ZBxi- 
ZB'x2 -M-F line (see the dotted line in FigO (b)). The F-M path runs parallel to the gold 
wires in the surface, crossing the 5x2 BZ through three different regions. The M-F line is 
perpendicular to gold wires. The surface/bulk and main atomic character of the different 
bands is identified by means of a MuUiken population analysis js^. 

Although a DZ basis is usually sufficient to obtain a quite good description of the occupied 
electronic states and the relaxed geometries in silicon systems, the use of a more complete 
basis set is necessary to describe the unoccupied part of the band structure even at low 
energies. For this reason all the band structures shown in the paper are calculated using 
a DZP basis set and slabs containing three underlying silicon bilayers (even if the relaxed 
geometry is obtained from a calculation using a DZ basis and/or a thinner slab). 

Finally, the Scannin g Tu nneling Microscopy (STM) images are simulated using the theory 
of Tersoff and Hamann 15511. 



III. RESULTS 



In this section we present our results for the different models of the Si(lll)-(5x2)-Au 
surface. We first focus on the energetics, relaxed geometries, and the electronic band struc- 
tures. We then turn our attention to the effect of the different silicon adatom contents and 
the simulated STM images, which we only analyze in detail for the most stable structural 
models. A summary of the relative energies of the calculated configurations, accompanied 
with a brief description of each of them, can be found in Table lU 

Before starting with the description of the results, it is interesting to point out some brief 
comments about the concentration of silicon adatoms on the Si(lll)-(5x2)-Au surface. A 
detailed study of the equilibrium situation has recently been performed by Kirakosian et 



al. |36, using STM. Their results indicate that, at equilibrium, only one adatom site 
is occupied out of every four possible sites, corresponding to a 5x8 adatom periodicity 
(if all the adatom sites were occupied we would recover a perfect 5x2 periodicity). The 
analysis of the adatom-adatom correlation functions obtained from the STM images reveals 
a strong suppression of those configurations with small adatom-adatom distances, a clear 
maximum corresponding 5x4 periodicities, and a long range oscillatory tail This was 
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interpreted in terms of a short range repulsion between adatoms plus a long range interaction 



term. In Ref. 



36l Kirakosian and collaborators showed that the density of adatoms can 



be increased by depositing additional amounts of silicon reaching an almost perfect 5x4 
arrangement of the silicon adatoms. Further deposition of silicon does not create a stable 
5x2 adatom structure. Instead the extra silicon atoms decorate the step-edge of the terraces 
on the surface. These observations seem to have at least two implications: (i) the optimal 
adatom concentration must be certainly lower than one adatom per 5x2 cell and, {ii) the 
structure of the reconstruction must be stable against relatively large changes of the content 
of adatoms ^62] since the density of silicon adatoms can be increased by a factor of two 
without, at least apparently, dramatic structural changes 

A systematic study of the energetics of the surface as a function of the adatom concentra- 
tion by means of first-principles electronic structure calculations is quite complicated. This 
is for two main reasons. First, the energies involved are rather small, which implies the need 
of very well converged calculations. A more serious limitation, however, is the necessity 
to use large supercells consistent with the low adatom densities. For this reason we have 
concentrated most efforts in the two limiting cases, involving respectively and 1 adatoms 
per 5x2 cell. The intermediate concentrations usually require drastic approximations. For 
example, Erwin [33] assumed that the main effect of the adatoms in the Si(lll)-(5x2)-Au 
surface is to dope the gold chains with electrons and studied the energetics of the system 
as a function of this doping. Here we go a step beyond and present explicit calculations 
for adatom contents down to 1/4 adatoms per 5x2 cell, consistent with a 5x8 periodicity, 
which indeed can be reached in experimental conditions j^^. Due to the large size of these 
systems we limit this study to our two most stable models, and only use the smaller DZ 
basis set. 



A. Marks and Plass model 

We start our investigation of the structure of the Si(lll)-(5x2)-Au surface using the 
model proposed by Marks and Plass [3^ from experimental data obtained with heavy-atom 
holography and high resolution electron microscopy. We use the label MP+ for this structure 
(see Table H}. The + superscript indicates that the structure contains silicon adatoms 
saturating some of the silicon dangling-bonds in the structure, what we call "conventional" 
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32, can be found 



silicon adatoms. A schematic view of this structure, as proposed in Ref. 
in Fig. ^ It has to be taken into account that, due to the hmitations of the experimental 
techniques, there are several uncertainties and assumptions in this structure. Only the 
atomic coordinates within the surface plane are accurate. The heights of the atoms over 
the substrate are only approximately resolved. The experimental beam error in combination 
with the size and complexity of the structure also limits the sensitivity to possible subsurface 
relaxations. As a consequence, the experimentally proposed structure only considers the 
reconstruction of the outermost bilayer and contains limited information about the registry 
between this surface bilayer and the underlying material. It is necessary to eliminate these 
uncertainties before one can undertake any serious study of the electronic structure of the 
MP"*" model. In order to do this while preserving all the information originally present 
in the MP"*" proposal, we started our study by performing constrained relaxations of the 
structure. The structure in Fig. ^ was relaxed using following degrees of freedom: (i) the 
height of the different layers and, (u) the lateral position of the surface layers with respect 
to the underlying bulk slab. The grouping of the atoms in different layers given in Ref. 
only implies approximately equal z-coordinates ( the z-axis is taken here along the surface 
normal). For this reason, in a second step, the atoms were allowed to relax in the z-direction 
while keeping their coordinates within the xy-plane. The resulting geometry preserves the 
bonding pattern of the original MP+ proposal, and provides a reasonable initial guess to 
start our search of the most stable models by performing full structural optimizations. 

We now consider in detail some of the structural patterns appearing in the MP"*" model of 
the surface. For this analysis we find useful the comparison with the Si(557)-Au surface, a 



Close 



related reconstruction that has been quite well characterized during recent years 



10|,|l2]- The stepped Si(557)-Au is formed after the deposition of ~0.2 monolayers 
of gold on vicinal (111), with the misorientation chosen along the [112] direction. The size 
and orientation of the terraces of the Si(557)-Au represent an analogous to the fiat 5x2 
unit cell but including a single silicon step ^3]- With half the gold coverage than the 
Si(lll)-(5x2)-Au surface, the terraces of the Si(557)-Au contain only one Au wire running 
parallel to the step edges. Gold atoms occupy silicon substitutional positions in the surface 
layer. This is supported both, by recent X-ray diffraction data and d_ensity functional 



calculations performed using a methodology similar to the one utilized here 



, which provide 



a consistent structural model of the surface. In particular, the highest stability of the 
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silicon substitutional sites for gold has been unambiguously demonstrated by the ab initio 
calculations. For example, the substitutional site was determined to be at least 1 eV/Au 
more stable that adatom-like positions, where gold sits on the surface saturating one of the 
silicon dangling-bonds, or even ~0.5 eV/Au more favorable than the adsorption decorating 
the step edges j8|. It seems, therefore, that the Au atoms on the Si(557)-Au surface exhibit a 
strong tendency towards three-fold silicon coordination. Gold atoms adapt to this situation 
without much strain, with typical Si-Au distances only a few percents larger than the bulk 
silicon bond length. 

In the light of these observations the bonding pattern of some of the gold atoms in the 
MP"^ model (Fig. seems quite peculiar. In particular, the gold atoms in the chain situ- 
ated at the left side of the "gold trench" (marked with an L in Fig. present a fourfold 
coordination. They are connected to three silicon atoms within the surface layer and, addi- 
tionally, to the silicon atom immediately below. Furthermore, the Si atoms neighboring to 
the mentioned gold atoms (see atoms a and a' in Fig. ^ present an unsaturated dangling 
bond which might be avoided with a slight structural change. 

It is interesting to note that the tendency of the gold atoms to occupy silicon substitu- 
tional positions in the top most layer cannot help to completely rationalize the structure. 
A three-fold bonding pattern of the gold atoms is inherently frustrated by the presence of 
a surface dislocation. In the MP"^ model this dislocation is located at the position of the 
right-hand gold wire (marked with R in Fig. P). Due to the change of the bonding sequence 
there are not three unpaired silicon electrons available for each of these Au atoms, but rather 
two. Therefore, they do not occupy a normal three-fold position and are quite likely to be 
displaced from the initial symmetric positions after relaxation. It should be noted that, 
in principle, the surface dislocation can be moved to different locations. In fact, we will 
see below that this provides a simple route to generate alternative structural models of the 
surface. 

The comparison between the structure of the Si(557)-Au reconstruction j^, 12 1 and the 



MP"*" model of Fig^ raises another interesting point. In the case of the Si(557)-Au surface 
the silicon atoms in the proximity of the step edge suffer a considerable rebonding. They 
form characteristic silicon structure which has been identified with the so-called "honeycomb 
chain" (HC) by several authors 0,0]. The presence of the silicon HC seems instrumental 
to understand the stability of the Si(557)-Au and related reconstructions (see, for example, 
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Ref. ISl). The silicon HC was initially proposed by Erwin and Weitering[38| as the main 
building-block of the (3x1) reconstruction induced on Si(lll) by the deposition of metals 
like Ag, Li, Na, K, Mg or Ba. The HC structure represents a large disturbance from the 
usual bonding pattern of silicon. The stability of the silicon HC stems fundamentally from 
the formation of a double-bond between two three-fold coordinated silicon atoms on the 
surface In the case of the (3x1) reconstruction a further stabilization mechanism comes 
from the fact that, after the transfer of the valence electrons from the metal atoms, the HC 
structure becomes electronically closed-shell. [6^ It seems somewhat surprising that the 



silicon HC structure, common to the (3x1) and Si(557)-Au metal induced reconstructions, 
is absent from the MP"*" model of the Si(lll)-Au-(5x2) surface. Indeed the MP"*" model 
seems to be based on an almost unreconstructed Si(lll) surface with a row of adatoms 
on top, and the more clear disturbance from this bonding pattern being the presence of a 
surface dislocation. 

We now proceed further with the structural relaxations of the MP"*" system. It is instruc- 
tive to focus first on a optimization were only the silicon degrees of freedom are taken into 
account. The gold atoms are constrained to remain at their initial coordinates. Due to the 
more directional bonding of silicon we can expect the structural changes to be simpler to 
analyze and somewhat less dependent on the particular choice of the initial guess in this 
case. Besides, as a stronger scatterer, we can assume that the gold positions to be better 
resolved in the experiment. The resulting geometry is plotted in Fig. IHfa). We observe two 
main effects. On the one hand, the HC configuration clearly emerges. One of the driving 
forces behind this rebonding is the movement forward of a and a' atoms in order to form 
an additional covalent bond with the silicon atoms in the underlying layer. The double- 
bonded "dimers" of the HC structure are formed by atoms h and c. This questions the 
location of the adatoms in the surface since, in principle, the dangling-bond associated with 
atom h' could disappear with the formation of a silicon double bond. We can observe the 
disturbing effect of the adatom on the HC structure. The appearance of the HC bonding 
pattern during the relaxation of the MP"*" structure confirms the results of recent density 
functional calculations by Kang and Lee who also made a geometrical optimization of 
the MP"*" model. The electronic bands calculated for this structure (not shown here) are 



also in quite good agreement with those presented by these authors in Ref. 



33. FigEl(a) 



also shows clearly what could be classified as a "stacking fault" in the structure (bonds of 
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atoms d and e coincide with those in the underlying sihcon layer). This stacking fault, which 
probably is energetically unfavorable, can be easily avoided by moving the position of the 
surface dislocation from the right-hand to the left-hand of the gold trench. Alternatively we 
can visualize this change (at least approximately) as a 180" rotation of the surface layer with 
respect to the underlying silicon structure. This transformation gives one of the structures 
discussed in the next section, which incidentally is almost identical to the "5x1" structure 



proposed recently by Erwin in Ref. 



I — 
33. 



When the relaxation of the MP+ system is continued without any constraints, the 
monatomic gold wires are strongly distorted as can be seen in Fig. |S1 (b). This distortion 
was not observed in the density functional calculations of Kang and Lee 37|. The reasons 
for this discrepancy are not completely clear at the moment. The break of the monatomic 
gold wires seems to be related with the presence of adatoms. If they are eliminated from 
the structure the gold atoms remain in two well separated parallel wires. Additionally, the 
strain introduced by the adatoms in the structure, results in the weakening of some Si-Si 
bonds in the surface layer (see the increased distance between atoms / and g). In spite of 
these strong structural distortions, the presence of adatoms in the structure is still energet- 
ically favorable as can be seen in Table HI These structural distortions are reflected in the 
band structure: following the nomenclature used in Ref. |37|, the band Si is shifted to higher 
energies respect to the S2 and a gap of ~0.3 eV is opened respect to the constrained case. 

In summary, our results suggest that neither the silicon structure nor the positions of the 
gold atom in the structure proposed by Marks and Plass are stable. Furthermore, in 
agreement with the general conclusions of Ref. 0, neither the STM images nor the band 
structure of the fully relaxed or the constrained relaxed MP"*" model seem to be in agreement 
with the experimental information. 



B. The Erwin models 



As discussed in the previous section, the MP model of the Si(lll)-(5x2)-Au surface recon- 
struction is characterized by the presence of a surface dislocation between one of gold wires 
and the neighboring silicon atoms. Other locations are possible for the surface dislocation. 
In particular, it can be translated to the other gold wire, this can also be assimilated to a ro- 
tation of the surface bilayer with respect to the underlying bulk silicon. This eliminates the 
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"stacking fault" commented in the previous section, and produces a new structural model. 
This structure is very similar to the "5x1" model recently proposed by Erwinjs^, and we 
refer to it as E(5xl). In Table|l]we can see that the E(5xl) model without silicon adatoms 
is slightly more stable than the relaxed MP structure. Fig. E](a) shows the relaxed structure 
of the E(5xl) model. The left (L) gold wire, where the surface dislocation is located, suffers 
a considerable dimerization, which is much smaller for the right (R) wire. The alternating 
Au-Au distances as obtained with VASP are, 4.06 A and 3.59 A for the L wire, and 3.82 A 
and 3.83 A for the R wire. The geometries obtained with SIESTA are very similar, specially 
those obtained with the more complete DZP basis set. Hereafter we name "SiAu complex" 
the structure formed by the two gold wires and the central silicon atom connecting them. 
The silicon structure in between two of such SiAu complexes is quite flat and resembles 
what could be described as a double honeycomb chain (DHC) silicon structure js^. The 
band structure along the direction parallel to the gold wires is shown in Fig. 0] (b). It shows 
several surface bands and has a metallic character. Those surface bands mainly associated 
with the Si-Au complex has been highlighted using solid symbols. Most of these bands are 
occupied and appear in the gap region. The unoccupied surface bands appearing in the 
gap are mainly associated with the silicon DHC. The most prominent feature is a dispersive 
band associated with the weakly dimerized (right) gold wire and the central silicon atom in 
the SiAu complex. This band is, in principle, metallic and close to half occupied. Although 
small gaps are opened associated with the crossings with other bands and slight geometrical 
distortions, it can be easily followed in Fig. |3](b) extending from ~1.3 eV below to ~2.3 eV 
above Ep. A similar band, with a similar origin, also dominates the band structure of the 



Si(557)-Au surface Q^yi- This band comes mainly from the sp^ lobes of the central silicon 
atom in the SiAu complex. There is also a strong hybridization with the 6p states of the 
gold atoms in the R wire. For this reason, they are better assigned to the Si-Au bonds 
connecting the central silicon with the R gold wire. Its large dispersion is due to the large 
overlap between these Si-Au bonds along the wire. The metallicity stems from the inability 
of gold (each gold atom only provides one valence electron) to saturate the bonds with all its 
silicon neighbors ^ . The other states in the Si-Au complex give rise to relatively fiat surface 
bands associated either with weakly overlapping silicon states or with the gold dimers. 



In Ref . 



33l it was also proposed that, under certain conditions, it could be energetically 



favorable to remove some of the over coordinated silicon atoms in the neighborhood of the 
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surface dislocation. Our relaxed structure for this model (hereafter E(5x2)) is shown in 
Fig. In this case both gold wires present an appreciable dimerization with alternating 
Au-Au distances of 4.37 A and 3.35 A for the left gold wire, and 4.16 A and 3.49 A for the 
right wire. Our SIESTA calculations with the smaller DZ basis set predict the E(5x2) model 
to be more stable, by at least 3.4 meV/A^, than both the E(5xl) model and the different 
variants of the MP model (see Table |ll). However, the difference between the E(5xl) and 
E(5x2) models is reduced with the use of more complete basis set. In particular, our plane- 
wave calculations predict both models to be degenerate within 0.1 meV/A^ (the E(5xl) 
slightly more stable). This agrees with the results of Ref. |33| where the E(5xl) model is 
predicted to be more stable than the E(5x2) variant by less than 1 meV/A^, and only after 
the addition of silicon adatoms the E(5x2) structure becomes favorable. 

The band structure of E(5x2) with zero adatom coverage is plotted in Fig. IHl The 
band st..uotu..e along tKe wi.es . in good ag,.een,ent with that reported in Ref. y for tMs 
structure. Again, the surface bands close to the Fermi energy come mainly from the SiAu 
complex. Like in the case of the E(5xl) model, the band structure is metallic. This is in 
disagreement with one of the latest and more detailed ARPES experiments which sug gest s 
that the Si(lll)-Au-(5x2) surface is a semiconductor with a band gap of at least 0.2 eV |29[. 
However, the metallic versus semiconducting character of this surface is still a matter of 
controversy. For example, the recent ARPES study by Himpsel and collaborators finds 
several metallic bands j30|. In fact, this reference and the scanning tunneling spectroscopy 
(STS) data of Ref. |3J indicate that the surface could be composed of alternate metallic 



and semiconducting regions along the gold wires. Our calculated band structure for the 
E(5x2) model is very close to being semiconducting. Just by shifting the Si band to higher 
energies by a few tenths of eV we could obtain a semiconducting surface. This might 
indicate that the metallic behavior is simply related to the limitations inherent to the local 
density approximation used here and the very simplified assumption that the monoelectronic 
eigenvalues can be directly identified with the photoemission peaks. In spite of its metallicity, 
several characteristics of the photoemission spectra are recovered by the band structure in 
Fig. ini The most prominent band observed experimentally starts at the boundary of the 
5x2 zone (ZBx2) dispersing downwards until it reaches a minimum at the boundary of the 
5x1 zone (ZBxi) js?, 28, 29, 3^. This band appears at binding energies between ~0.2 eV 
and ~1.3 eV. Following Erwin j33|, we can try to identify this band with our S2 band, whose 
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maximum appears close to Ei? in the neighborhood of ZBx2- However, it becomes difficult 
to follow the dispersion of this surface band as we move to higher binding energies for two 
reasons: i) the band enters the region of the projected bulk bands, becoming a surface 
resonance and, ii) other surface bands coming from the same region of the surface appear 
in the energy interval between -0.5 and -1.2 eV. This last point is widely consistent with 
the experimental data in Ref. |29|, where three additional bands are identified for binding 
energies larger than 0.5 eV. 

Losio and collaborators [3] reported an interesting effect, a continuous dimensionality 
transition of the main surface band. The character changes from strongly one-dimensional 
at the band maximum (i.e. only dispersing in the direction parallel to the gold wires) to two- 
dimensional at its minimum (i.e. with a non- negligible dispersion also in the perpendicular 
direction). The strong one-dimensional character of the surface states close to Ei? has also 
been eonfir.ed the ,„ost .eeent ARPES ,„easu_ QQ. The dispersions in the 
direction perpendicular to the wires can be found in Fig. IHl (b) and (c). The band widths 
are rather small for most surface bands. An effect similar to the reported dimensionality 
transition can be seen in the case of the Si band. It is tempting to assign the experimentally 
observed effect to the S2 band (see the different dispersion of bands S2 in panel (b) and S' 
in panel (c)). However, as commented above it is not so simple to follow the S2 band as it 
disperses downwards. In fact, we can locate what seems to be an avoided crossing between 
the S2 and the S3 bands half way along the ZBx2-ZBxi path in Fig. IHl (a). Therefore, we 
think that the S' band in panel (c) is rather related to the S3 band than to the S2 band, 
and the dimensionality change would be absent from our results. Also the energy position 
of the band S' (~-0.5 eV) is quite far from the ~-1.3 eV found experimentally for the 
band minimum. Therefore, in contrast to Erwin 3^ we conclude that our calculated band 
structure for the E(5x2) model does not provide a direct explanation to the observation by 
Losio et ai. 

Similarly to the surface bands of the E(5x 1) model, the Si, S2 bands in Fig.|ni(a) have the 
largest weight in the central Si atom in the SiAu complex. The Si band can be associated 
with SiAu bonds connecting the central Si with the left gold wire. This SiAu bonds have a 
small overlap and this is translated in a quite fiat band. The two dispersive S2 and S3 bands 
have a stronger weights in the other SiAu bonds, which have a larger overlap and, therefore, 
present a stronger dispersion. 
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We now explore the role of the silicon adatoms in these structures. We first studied 
the stability of the adatoms in the E(5xl) model when they are located over the silicon 
part of the surface reconstruction, i.e on sites equivalent to those occupied by the adatoms 
in the original MP proposal. It is interesting to note that the role of the silicon adatoms 
in such positions is indeed not very clear. The stability of the adatoms in typical silicon 
reconstructions stems from the fact that each adatom can saturate three dangling bonds 
in the surface at the expense of creating just an additional dangling bond. The energy 
gained in this process usually overcomes the strain energy caused by the addition of the 
adatoms. However, the E(5xl) model in Fig. |3] (a) does not have silicon dangling bonds. 
The appearance of unsaturated dangling bonds is avoided by the formation of the double- 
bonded silicon dimers that characterize the HC configuration. In fact, the only metallic 
band in this model comes from the SiAu complex as explained above. In accordance with 
these observations, we found extremely difficult to reach a stable configuration, i.e. with all 
the components of the forces below our threshold, for the silicon adatoms over the silicon 
sites of the E(5xl) model. Finally, after several hundreds of optimization steps this model 
spontaneously relaxed into a new structure. This structure, labeled N"*" in Table. HI belongs 
to a new family of structural models for the Si(lll)-Au-(5x2) surface found in this work for 
the first time and described in more detailed in the next subsection. 

In the light of the previous comments, a more stable adsorption site for the silicon adatoms 
would be on top of the SiAu complex. This has been previously proposed by Erwin j^^, 
and is confirmed by our calculations. Table H] shows the changes in surface energy after the 
addition of one silicon adatom per 5x2 unit cell. The behavior is opposite for the E(5xl) 
and E(5x2) models, with the addition being energetically favorable for the later model. 
The E(5xl) remains metallic after the addition of the adatom, and the dispersive band 
associated with the SiAu complex remains quite unchanged. The situation with the E(5x2) 
model is different. In agreement with the results in Ref. I^^we find that the band structure 
becomes semiconducting after the addition of the adatoms. The corresponding atomic and 
electronic structure can be found in Fig. [3 (a) and (b) respectively. The surfaces bands with 
a larger contribution from the atoms in the SiAu complex has been highlighted using solid 
symbols. It has been impossible to identify a band that can be solely assigned to the adatoms. 
We can see that the band structure of the E(5x2) suffers major modifications after the 
addition of adatoms, at least for the large concentrations considered here. Besides the fact 
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that the structure becomes semiconducting, the agreement with the detailed photoemission 
experiments of references |29| and |30| seems to be somewhat degraded. 



C. New structural model 

In this section we present a novel structural model for the Si(lll)-(5x2)-Au surface 
reconstruction that has been found during our investigation. Our slab spontaneously relaxed 
to this new structure while trying to optimize a modified version of the E(5xl) model 
commented in the previous section. The new structure can be found in Fig. |H1 and will be 
referred here as model N. Table |l] shows that the energy of the new model compares favorably 
with those of the other structures proposed to date. In fact, within our calculational scheme 
it is the most favorable structure. The difference with the second most stable model without 
adatoms, the E(5x2), is 4.7 meV/A^. This difference is reduced to 4.1 meV/A^ when using a 
more complete DZP basis set as shown in Table ITTl These energy differences are quite small, 
so further studies have been performed in order to drive more definitive conclusions. First, we 
have repeated our calculations using the PBE jsTj GGA exchange and correlation functional 
instead of LDA. The new structure continues to be more stable by 3.6 and 5.1 meV/A^ using, 
respectively, a DZ and a DZP basis set. As a second step, the energy ordering between the 
N and the E(5x2) structures has been confirmed using VASP and slabs containing three 
and four silicon double-layers. The new model is more stable than the E(5x2) by at least 
2.6 meV/A^. These results convincingly establish, at least within the framework of density 
functional calculations, the larger stability of our new structural model compared to previous 
proposals in the limit of negligible adatom coverage 

Given the small energy differences between both models it may be interesting to esti- 
mate the effect of the vibrational degrees of freedom in the surface free energy 7(T). The 
vibrational contribution can affect the energy ordering even at zero temperature due to 
the zero-point energy, and its importance grows with temperature. Unfortunately, an ac- 
curate estimation of the vibrational surface free energy 7«b(T) is a formidable task that 
would require the detailed calculation of the dynamical properties (phonon band struc- 
ture) of the different surface models. This is a computationally very demanding calcu- 
lation that is beyond the scope of the present paper. We can obtain a rough estima- 
tion of the vibrational contribution to the difference of the surface free energies between 
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the different structures A7(T) following Ref. |58l We have A7(T) = AEgurf + A'jmbiT), 



where A7„.(T) ^ 3N§''''\F{T,ujfSf^) - F{T,u;,uik)] - 3Ng[F{T,u;Z^ - F(T,^,„,fc)] + 



3Nau[F{T,uj^I^^'^^) — F{T,uj^^)]. Here AEgurf is given in Table HI and is independent of 
the temperature T; Ng^^^'^^ and Ng^ are the number of silicon atoms per unit cell in both 
surface structures, and A^^^ the number of gold atoms; F(T, cu), given in the Appendix of 
Ref. |58|, is the free energy of a given vibrational mode u; the frequencies oo^urfi ^surf'^\ 
and ujbuik characterize the average vibrational properties of the silicon atoms in both sur- 
face structures and in bulk silicon, while cu^^^^^^ and uj^^ those of the gold atoms in both 
surfaces. We take for Ubuik values in the range of 50-70 meV, and Ugurf ranging from 0.5 to 
1.5 the ujbuik value. Within these range of parameters, if cuf^^j and oof^^j'^^ differ less than 
a ~10%, then A7^jfc(T) stays within ~ ±2 meV/A^ for temperatures up to 300 K. If the 
vibrational properties of both surface models differ more significantly, then A7^,jf,(T) can af- 
fect the relative order of the structures at much lower temperatures. However, we should not 
expect strong differences in the average vibrational frequencies of the E(5x2) and N models. 
Both models present very similar bonding patterns and structures. It is interesting to notice 
that A'~fmb{T) is nonzero even if the vibrational properties of both structures are identical, 
i.e. ujf}^^'^^=uj^^^f =ujsurf and 00^1^^"^^ =uj^^. This reflects the different number of silicon 
atoms in the unit cell of the two surface reconstructions. In this case we have A'-)mb{T) ^ 
— Ng-)[F{T,uJsurf) — F{T,uJbuik)]- Using the same parameters as above we ob- 
tain A'-fmb{T) within ±1.5 meV/A^ up to ~1000 K. Thus we can conclude that the energy 
ordering obtained in the present total energy calculations is not altered by the vibrational 
contribution to the free energy up to, at least, room temperature. 

In the new structure the gold wires along the [110] direction present a dimerization 
comparable to the E(5x2) structure. The alternating Au-Au distances are 3.24 A and 
4.40 A (3.19 A and 4.45 A) along the right (left) wires. The distance between nearest 
neighbor Au wires along the [112] direction is smaller in the N structure (3 A) than in 
the E(5x2) structure (3.8 A). The later value b eing in better agreement with the ~3.9 A 
deduced from the HREM studies of the surface. 

Similarly to the E(5xl) structure, most of the surface of the N model is covered with a 
silicon double honeycomb chain structure jssi] . One of the silicon atoms in the DHC appears 
at a higher position over the surface. This indicates that this atom has a charged dangling- 
bond and, therefore, is trying to develop a sp^ hybridization. This atom is expected to be 



20 



more visible in the STM images and to provide a preferential site for adsorption on the 
surface, in particular for possible silicon adatoms. The boundaries between the DHC stripes 
are occupied by the SiAu complex, in which a central silicon atom appears bonded with 
three gold dimers. 

The band structure of the new structure is plotted in Fig. El The general features are in 
good agreement with the most recent ARPES studies 0,01; although some of the details 
are different. The most dispersive and prominent surface bands are quite similar to those 
found for the E(5x2) model. The surface is predicted to be semiconducting, which agrees 
with the results of Ref. 13 The bands named Si and S2 by Matsuda et al. j^l can be easily 
identified in our calculation, and we use the same notation. Other less dispersive surface 
bands are also observed in our calculated band structure. These can be tentatively identified 
with those labeled S3 and S4 by Matsuda et al. However, the S3 band appears shifted to 
lower binding energies by a few tenths of eV. We can relate this upward energy shift to 
the use of the LDA in our calculations, which is likely to be less suited to describe more 
localized (less dispersive) states. Besides this energy shift, the sole major discrepancy with 
the experimental band structure in Ref. is the absence of the S3 band. However, this 
band is not so clearly resolved in the experiments as the others. 

Different symbols are used in Fig. El according to the main atomic character of the bands. 
Si and S2 come from the Si-Au bonds in the surface (solid symbols). This is common to 
most of the models studied in this paper: the most dispersive surface bands always originate 
in the Si-Au bonds, with the main character corresponding to the 3p states of the central 
Si atom, and a strong hybridization with the Qp states of the neighboring Au atoms. The 
flat S3 band corresponds to the silicon dangling bonds in the middle of the DHC structure 
(open triangles). The also quite fiat S4 band is mainly associated with the bonds between 
the gold atoms and the silicon atoms in the border of the DHC structure (open squares). 
We find several unoccupied surface bands whose atomic character is difficult to determine. 
One of these bands is located at energies very close to Ei?, particularly near the F point. 
The metallic/semiconducting character of the surface is thus governed by the position of 
this band. This situation is very similar to that already observed for the E(5x2) model, 
although in this case the band reaches to lower energies and becomes partially occupied 
driving the system to metallic. 

In agreement with experiment, most surface bands show a strong ID character in our 
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new structural model as can be seen in Fig. 121 (b) and (c). This is particularly clear in panel 

(b) , where most states are located within the bulk gap. In the region displayed in panel 

(c) (at the zone boundary of the 5x1 Brillouin zone) the S2 and Si bands merge with the 
bulk bands, becoming surface resonances. It is no longer possible to identify the Si and 
S2 resonances with a single band of our finite slab and, as a consequence, it is difficult to 
follow the band dispersion of these spectral features in the direction perpendicular to the 
gold wires. However, from the data in panel (c) it is clear that the combined effect of the 
possible dispersion, plus the broadening of the resonances extends over a range of ~0.2 eV, 
much larger than its dispersion for ener gies closer to Ep. This is broadly consistent with 
the ID to 2D transition reported in Ref. |28| for the most prominent photoemission feature 
as the binding energy increases. 

We now explore the structure and energetics of the model N under the addition of one 
silicon adatom per 5x2 unit cell. We tried several different adsorption sites: directly on the 
SiAu wire following the proposal by Erwin [2^ (referred as N*), and bonded to the prominent 
dangling bond in the DHC structure occupying hollow H3 (N"*") or top T4 (N"*"') sites j^. 
As shown is Table HI this high coverage of adatoms is energetically unfavorable in all cases 
by at least 2.7 meV/A^. This is in contrast with the situation for the E(5x2) model, where 
the addition of one silicon adatom per unit cell is slightly favorable. In the N* structure (not 
shown) the silicon adatoms tend to locate in a peculiar bridge position between two gold 
dimers along the [110] direction. The structure of the N"*" model is shown in Fig.^l(a). The 
silicon atoms bonded to the adatom adopt a typical silicon configuration although, contrary 
to what is observed for the clean Si(lll) surface, the hollow site is preferred over the top 
site Q. 

The band structure of the N"*" surface is shown in Fig. (b). It is very similar to that 
found for the model without adatoms. The Si and S2 are largely unchanged, which clearly 
indicates its origin in the SiAu complex. The fiat S3 band disappears from the gap region as 
a consequence of the saturation of the dangling bond with the adatom. A new unoccupied 
band, associated with the adatoms, appears instead. This new band can be found around 
~0.6 eV above in Fig. HU] (b). 
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D. Adatom coverage 



So far we have only considered the hmiting cases with zero or maximum adatom coverage, 
which correspond to a number x of sihcon adatoms per 5x2 unit cell equal, respectively, to 
and 1. However, the experimental evidence indicates that the equilibrium concentration 
is X ~ 1/4, corresponding to one adatom per 5x8 supercell. Under silicon rich conditions 
the adatom coverage can be increased in the experiment only up to a; ~ 1/2, consistent 
with a 5x4 periodicity. We have performed explicit calculations for x = 1/2, x = 1/3, and 
X = 1/4 for our two most stable models of the reconstruction in order to simulate these 
situations that can be reached experimentally. Due to the very large supercells necessary for 
these calculations (up to 273 atoms), we have performed them with the SIESTA code and 
restricted to the use of a DZ basis set for silicon. The results of the energetics as a function 
of the adatom content can be found in Table UTTl and in Fig. ^2 The behavior is opposite for 
both models, N and E(5x2). It should be kept in mind that model N favors the adatoms in 
hollow sites over the silicon surface, while in the E(5x2) structure the adatoms sit on the 
gold chains are more favorable. 

The surface energy monotonously decreases as a function of the number of adatoms for 
the E(5x2) model. We do not find any evidence of an energy minimum as a function of the 
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adatom concentration. This is in contrast with the suggestion made by Erwin in Ref. 
that reference the addition of adatoms was studied using the following simplification: it was 
assumed that the sole effect of the adatoms is to dope the gold chains with electrons and the 
energy of the system was studied as a function of the doping. Erwin found a minimum of the 
total energy for 0.5 extra electrons per 5x2 unit cell. Since each adatom was found to donate 
two electrons to the surface, this would correspond to the observed adatom concentration 
at equilibrium of x ~ 1/4. However, our simulations introducing explicitly the adatoms 
in the structure do not confirm this behavior. The surface energy of the E(5x2) structure 
always decreases as the adatom concentration is increased. However, the slope of the curve 
becomes very small for intermediate adatom concentrations, showing a weak dependence of 
the surface energy in that region. We cannot completely rule out the presence of a minimum 
for the surface energy at very low adatom concentrations. However, it seems quite unlikely 
looking at Fig ^2 It could also be argued that the DZ basis set is not flexible enough to 
produce the correct behavior. This seems quite improbable looking at the data in Table |h1 
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which clearly show that the energy changes induced by the addition of adatoms are weakly 
dependent on the details of the calculation. 

In the case of the new model N, the Fig. II II shows that the surface energy systematically 
increases as a function of the adatom concentration. With the DZ basis set the N model 
is always more stable than the E(5x2) structure. Using a more complete basis set and a 
converged plane-wave calculation we find a crossing: the new model is always more stable at 
low adatom coverage, but becomes unstable compare with the Erwin "5x2" model at larger 
coverages. Scaling the data calculated with the DZ basis set to reproduce the VASP results 
at the end points (i.e. x = and x = 1) we can estimate that the crossing occurs at x ~l/2. 
We can conclude then that the N model is, at least in the framework of density functional 
calculations, more favorable than the E(5x2) model for adatom concentrations below ~l/2 
adatoms per 5x2 cell. 



E. Simulated STM images 

The STM images of the Si(lll)-(5x2)-Au surface are characterized by the presence of 
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underlying lattice 

to silicon adatoms However, the origin of the "Y" -shaped features is less clear. 

Figures IT^ and IT!^ present our simulations of the STM images for the E(5x2) model at 
-0.8 eV sample bias and the N model at -0.6 eV, respectively. The simulations have been 
performed for a 5x4 arrangement of the silicon adatoms, corresponding to concentration of 
adatoms that can be actually reached in the experiment. In agreement with Ref. |33| and 
the experiments the silicon adatoms show as very pronounced bright protrusions. With 
the adatoms directly sited on the gold chains, the bright spots appear in the middle of the 
underlying row structures for the E(5x2) model. For the N model they appear in a more 
lateral position. This seems to be in somewhat better agreement with the experimental 
images (see, for example, the Figure 1 (b) in Ref. |56|). "Y"-shaped features can be identified 
in the simulated STM images of both N and E(5x2) models. The possible candidates have 
been highlighted in the Figures IT2l and IT3l (see also Ref . l33l) . The identification is, however, 
more clear in the case of the less symmetric E(5x2) structure. 
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IV. CONCLUSIONS 



We have performed a systematic study of different models of the Si(lll)-(5x2)-Au sur- 
face reconstruction by means of first-principles density-functional calculations using the 
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40| and the VASP 
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43 1 codes. We start our investigation with the struc- 

n 

tural model proposed by Marks and Plass [3^ (MP). This is the most detailed model of 
this surface reconstruction solely based on experimental information to date. Therefore, it 
provides a logical starting point for our study. We have also considered different variants of 
the relaxed MP model, including the structures recently proposed by Erwin jsS^, and a new 
structure found during our simulations. Within the computational schemes used here this 
new structure is the most favorable energetically, at least in the regime of low concentration 
of silicon adatoms. In general, we find a reasonable agreement between our results and 

energy differences 

between different models are quite small, with most structures lying in a narrow range of 
surface energies of less than 10 meV/A^ (the estimated error bar for our energies is of the 
order of 1-2 meV/A^). This, together with the uncertainties arising from the use of the 
local approximation to the density functional theory, make difficult to draw definitive con- 
clusions solely based on the energetics. The comparison of the calculated band structures 



and local density of states, respectively, with the available ARPES data 
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and the 



STM images j23|, |2J, |3J] becomes then instrumental in order to identify the most plausible 



candidates for the equilibrium structure. In the following we summarized some our main 
conclusions: 



i) Like in the case o: 
silicon surfaces la 



the reconstructions formed by the deposition of gold on stepped 



structure emerges as a 



the silicon honeycomb chain (HC) 
fundamental building block of the reconstruction. In agreement with the result of Hang and 
Lee js^] , the silicon HC is formed spontaneously during the relaxation of the MP model. The 
HC is also present in the optimized geometries of all the other structural models considered 
in our work. 

ii) For the MP model we agree with the main conclusions of Ref. that neither the 
simulated STM images nor the calculated band structure compare satisfactorily with the 
experimental data. 
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in) We have studied in detail the models proposed in Ref. by Erwin, the E(5xl) and 
E(5x2) structures. The E(5xl) model is quite similar to the MP structure: they correspond 
to two possible positions, at opposite sides of the SiAu complex, of the surface dislocation 
present in these structures. The E(5xl) model and its E(5x2) variant are energetically 
degenerate at zero adatom coverage. However, these two structures show a different behavior 
against the addition of silicon adatoms: it is always unfavorable for the E(5xl) model, while 
tends to increase the stability of the E(5x2) model. 

iv) We have explored a different position of the surface dislocation: at the center of the 
SiAu complex. We arrive in this way to a new structure, the N model. According to our 
calculations this new structure is more stable, at least for low coverages of silicon adatoms, 
than any of the models proposed to date. The distance between the gold wires in this model 
is ~3 A, which seems somewhat small compared to the ~3.9 A deduced from the HREM 
measurements js^. 

v) The calculated band structures of the E(5x2) and N models without adatoms are 

n 

aite similar and appear to be in reasonable agreement with the available ARPES data ^81 
n 
, QDI|. The other models fail to reproduce the main features observed experimentally. 

The agreement seems to be particularly good in the case of the N model. According to our 

analysis the most prominent and dispersive surface bands, named Si and S2 in Ref. |29|, come 

from the atoms in the SiAu complex. In the case of the N model the silicon adatoms tend 

to adsorb on the silicon part of the surface, i.e. bonded to three silicon atoms in the surface 

layer. As a consequence, the topology and the energy position of these bands are quite 

insensitive to the coverage of silicon adatoms. This contrast with the situation found for 

the E(5x2) model. Here the silicon adatoms tend to adsorb directly on the SiAu complex, 

thus causing a notable modification of the surface bands that worsens the agreement with 

the experimental ARPES spectra. 

vi) We have studied the energetics of the E(5x2) and the N models as a function of the 

n 

concentration of silicon adatoms. Contrary to the suggestion of Ref. we do not find 
any evidence of a minimum of the surface energy of the E(5x2) model as a function of 
the adatom coverage. The surface energy always decreases with the addition of adatoms, 
although the changes are very small in the range of x between 1/2 and 1/4, where x is the 
number of adatoms per 5x2 unit cell. For the N model the addition of the adatoms is always 
unfavorable. As a consequence of this opposite behavior, the E(5x2) structure becomes 
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more stable than the N structure in the hmit of relatively large adatom concentrations, 
for (notice that x = 1/2 corresponds to a 5x4 periodicity). According to this 

picture the exact content of adatoms is instrumental to determine the equilibrium structure 
of the reconstruction within the range of experimentally realizable adatom coverages. This 
introduces a new degree of complexity that should be taken into account when analyzing the 
experimental information. In particular, this might be behind the observed phase separation 
into 5x4 and 5x2 patches jsollsg. 

vii) The simulated STM images of the most stable models, N and E(5x2), are in broad 
agreement with the experimental images. The silicon atoms produce bright spots which are 
located in the middle of the underlying row structures for the E(5x2) and in a somewhat 
more lateral position for the N model. In both cases "Y" -shaped features similar to those 
observed in the experiment can be found. However, they are more clear in the case of the 
E(5x2) model 33| where the structure surrounding the gold chains is less symmetric. 
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Model 


Description 


AEsurf (meV/A2) 


MP+ 


Marks and Plass model after a constrained relaxation 


+46.8 


RMP+ 


Fully relaxed MP"*" structure 


+5.4 


RMP 


Relaxed MP"*" structure without adatoms 


+8.3 


E(5xl) 


Erwin "5x1" 


+4.8 


E(5xl)* 


E(5xl) with adatoms on the Au- wires 


+6.5 


E(5x2) 


Erwin "5x2" 


+1.4 


E(5x2)* 


E(5x2) with adatoms on the Au- wires 


0.0 


N 


New model 


-3.3 


N+ 


N with silicon adatoms in H3 positions 


-0.6 


N+' 


N with silicon adatoms in T4 positions 


+1.2 


N* 


N with adatoms on the Au-wires 


+2.4 



TABLE I: Summary of the structural models studied here for the Si(lll)-(5x2)-Au reconstruction 
and their relative surface energies (AEgurf)- Those structures containing adatoms have one silicon 
adatom per 5x2 cell, i.e. the concentration of adatoms is maximum. Superscript + indicates the 
presence of "conventional" adatoms saturating silicon dangling bonds in the surface. Labels H3 and 
T4 refer, respectively, to adatoms occupying hollow and top sites 5^. The presence of adatoms 
located on top of the Au wires is indicated by a * superscript. The data in this table have been 
calculated using the SIESTA code with a DZ basis for silicon and DZPs-SZd basis for gold. The 
slabs contained two silicon bilayers below the surface layer (see Fig.^(a)). All energies are referred 
to that of the structure recently proposed by Erwin in Ref. |33 . 
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Model 


AE,„,y (meV/A2) 




SIESTA 


VASP 




DZ-3 blys 


DZ-2 blys 


DZP-2 blys 


2 blys 


E(5x2) 


+1.3 


+1.4 


+1.6 


+1.4 


E(5x2)* 


0.0 


0.0 


0.0 


0.0 


N 


-3.4 


-3.3 


-2.5 


-1.3 


N+ 


-1.0 


-0.6 


+0.9 


+1.5 



TABLE II: Convergence of the relative surface energies (AE^ur/) of the most stable structural 
models respect to the basis set and the thickness of the slabs used in the calculations. The first 
column shows the data obtained with the SIESTA code using a DZ basis for silicon and three 
silicon bilayers below the surface to construct the slab. In the second column a slab with only two 
underlying silicon bilayers was used. The third and four columns are obtained using the thinnest 
slab and, respectively, a DZP basis set for silicon and the VASP plane-wave code. 



Model 




AE 


/ (meV/A2) 






X = 


X -g 


X — 3 


X — 2 


X = 1 


E(5x2)J 


+1.39 


+0.90 


+0.86 


+0.81 


0.0 


N+ 


-3.35 


-1.98 


-1.72 


-1.23 


-0.60 



TABLE III: Relative surface energies (AEsu^.j) of the most stable structural models as a function 
of X, the number of silicon adatoms per 5x2 unit cell. The calculational parameters here are the 
same used in TableUl Notice that E(5x2)*.^o=E(5x2) and N+^o=N. 



FIG. 1: (color online). Schematic view of a typical slab used in our calculations. It shows the 
model proposed by Marks and Plass (MP"*") [3j| for the Si(lll)-(5x2)-Au surface reconstruction. 
Large circles in the surface layer represent the gold atoms. The bottom surface of the slab is 
saturated with hydrogen atoms, (a) Side view and (b) top view with some of the silicon atoms in 
the surface and the two gold chain labeled (see the text). 
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FIG. 2: (a) Comparison of the bidimensional Brillouin zones corresponding to (5x1), (5x2) and 
(5x4) supercells on the Si(lll) surface, (b) Schematic view of the path (dotted hnes) used to plot 
the band structures in this work (F- ZBx2-ZBxi-ZB'x2"M"r). Its relation with the (5x1) (dashed 
hnes) and (5x2) (sohd lines) Brillouin zones is indicated, and some special points are defined. 

FIG. 3: (color online), (a) MP+ after the optimization of the position of the silicon atoms in the 
structure. The gold atoms are kept in the positions obtained after the initial contrained relaxation 
of the experimental coordinates. The silicon honeycomb chain (HC) structure has been highlighted, 
(b) The same structure after full relaxation (RMF"*"). See the text for the labels of the different 
atoms. 

FIG. 4: (color online), (a) Relaxed geometry of the E(5x 1) model with zero adatom coverage and, 

(b) the corresponding band structure. Solid symbols indicate those bands with a larger weights in 
the atoms of the SiAu complex. The energies are referred to the Fermi level. 

FIG. 5: (color online). Relaxed geometry for the E(5x2) model. 

FIG. 6: Band structure of the E(5x2) model with zero adatom coverage parallel (r-ZBx2-ZBxi- 
ZBx2 path in (a)) and perpendicular to the gold wires i (F-M path in (a) and panels (b) and 

(c) ). Surface bands with the larger contributions coming from the atoms in the SiAu complex are 
indicated by filled circles. The energies are referred to the Fermi level. 

FIG. 7: (color online), (a) Relaxed geometry of the E(5x2)* model (containing one adatom per 
5x2 unit cell) and, (b) the corresponding band structure. Surface bands with a larger weight in 
the atoms of the SiAu complex are marked with filled circles. The energies are referred to the 
Fermi level. 

FIG. 8: (color online). New structural model for the Si(lll)-(5x2)-Au reconstruction. This is the 
most stable configuration of the surface according to our calculations (see Table DJ). 
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FIG. 9: Band structure corresponding to the N model with zero adatom coverage parallel (F- 
ZBx2-ZBxi-ZB'x2 path in (a)) and perpendicular to the gold wires (F-M path in (a) and panels (b) 
and (c)). Surface bands are marked according to its main atomic character: filled circles indicate a 
strong contribution from the atoms in the SiAu complex, open triangles from the silicon dangling 
bonds in the middle of the double honeycomb chain (DHC), and open cubes from those silicon 
atoms at the boundaries of the DHC stripes, neighboring to the gold wires. The energies are 
referred to the Fermi level. 

FIG. 10: (color online), (a) Optimized geometry of the model (containing one adatom per 
5x2 unit cell) and, (b) the corresponding band structure. Surface bands with strong contributions 
from the SiAu complex are indicated by filled circles, while the states associated to the adatom are 
marked with open triangles. i The energies are referred to the Fermi level. 

FIG. 11: Relative surface energies as a function of the adatom content. Explicit calculations 
have been performed for several adatom concentrations using the smaller DZ basis set (circles). 
The results obtained with the DZP basis set (diamonds) and with plane-wave VASP calculations 
(triangles) for the two limiting cases are also shown for comparison. All energies are referred to 
those of the E(5x2)* model. 

FIG. 12: (color online). Simulated STM image of the E(5x2) model with a sample bias of -0.8 eV 
and an adatom concentration corresponding to a 5x4 periodicity. A possible candidate for the 
Y-shaped structure is schematically indicated. The atomic structure is superimposed with the 
simulated image in the lower part of the figure. Large circles indicate the positions of Au atoms. 

FIG. 13: (color online). Same as Fig. ^Jbut for the new N structure and a sample bias of -0.6 eV. 
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